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INTRODUCTION 


This  report  describes  our  technical  accomplishments  during  the 
first  six  months  of  the  present  contract.  The  report  is  in  two  parts. 

In  Part  A  the  results  of  a  one-dimensional  progression  failure 
modeling  of  rock  static  behavior  are  described.  In  Part  B  our 
experimental  findings  in  the  mechanical  behavior  of  rock  under  compressive 
cyclic  fatigue  are  detailed.  Part  B  is  to  be  presented  at  the  13th 
Symposium  on  Rock  Mechanics,  Urbana,  Illinois,  September  1,  1971. 

Work  is  now  underway  in  rock  tensile  cyclic  fatigue.  A  method 
has  been  found  to  test  this  fatigue  behavior  and  preliminary  results 
suggest  that  the  concept  of  a  fatigue  limit  for  rock  in  tension  is  viable. 

A  full  account  of  our  progress  in  this  aspect  of  the  program  will  be 
part  of  our  next  technical  report. 
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A  ONE-DIMENSIONAL  PROGRESSIVE  FAILURE  MODEL  OF  ROCK 


Robert  G.  Lundquist 
Robert  W.  Heins 


ABSTRACT 

It  is  well  known  that  the  strength  and  deformation 
properties  of  rock  are  highly  variable.  It  has  also  been 
observed  that  rock  does  not  fail  in  a  weakest  link  manner  but 
progressively,  maintaining  considerable  load-bearing  capac¬ 
ity  beyond  its  peak  load.  It  is  hypothesized  that  the  vari¬ 
ability  of  properties  and  progressive  failure  behavior  are 
direct  consequences  of  the  statistical  inhomogeneity  of  tho 
material.  A  one-dimensional  model  is  used  to  demonstrate 
that  progressive  failure  can  Indeed  resu”  from  inhomogeneity 
of  the  material  and  to  show  what  the  effects  of  varying 
statistical  distributions  of  properties  might  be  on  the  gross 
deformation  and  failure  of  rock. 
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"You  see,  of  course,  if  you're  not  a  dunce, 
How  it  went  to  pieces  all  at  once,  — 

All  at  once,  and  nothing  first,  ~ 

Just  as  bubbles  do  when  they  burst. " 

The  Deacon's  Masterpiece 
Oliver  Wendell  Holmes 


CHAPTER  I 
INTRODUCTION 


Problem  Statement 

In  the  Interpretation  of  laboratory  tests  and  in  the  design  of 
structures  in  rock,  we  often  assume  that,  like  the  Deacon's  "wonderful 
one-hoss  shay,  "  rock  fails  "all  at  once,  and  nothing  first.  "  It  is  easily 
demonstrated  that  rock  does  not  fail  in  this  way  under  most  practical 
loading  conditions,  and  it  is  the  purpose  of  this  study  to  show  the  con¬ 
sequences  of  this  misleading  assumption  and  to  present  a  simple  model 
for  rock  failure  which  allows  the  interpretation  of  rock  strength  tests  and 
the  design  of  rock  structures  to  be  put  on  a  more  rational  basis. 

There  are  several  Implications  in  the  idea  of  failing  "all  at  once 
and  nothing  first"  which  we  must  examine.  The  first  and  most  obvious 
Is  that  once  the  strength  of  the  material  is  exceeded  by  the  stress,  it  is 
no  longer  capable  of  sustaining  load.  An  abundance  of  recent  evidence 
has  shown,  however,  that  in  compressive  and  mixed  stress  fields,  and 
in  perhaps  every  loading  except  uniaxial  tension,  most  rock  can  sustain 
considerable  fractions  of  its  peak  load  after  this  load  has  been  exceeded. 
This  is  the  concept  of  the  "complete  stress- strain  curve,  "  perhaps  more 
precisely  called  the  "complete  load-displacement  curve"  since  only  mean 
stress  and  strain  can  be  determined  after  failure.  It  is  this  behavior 
which  led  Jaeger  and  Cook  (1)  define  the  brittle  state  of  a  material  as 
that  condition  in  which  "it  ability  to  resist  load  decreases  with  increas¬ 
ing  deformation.  " 
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A  second  implication,  which  will  be  seen  to  be  intimately  con¬ 
nected  with  the  first,  is  the  idea  of  the  "weakest  link.  "  By  analogy  to 
the  ancient  adage  that  a  chain  is  no  stronger  than  its  weakest  link,  it  is 
assumed  that  if  failure  ensues  from  exceeding  the  strength  of  any  part  of 
a  rock  structure,  the  entire  structure  will  fail.  Failure  must  here  be 
defined  as  the  loss  of  the  entire  load  resisting  ability  of  the  structure. 

To  see  that  this  is  untrue,  we  can  examine  some  physical  observations 
which  have  accompanied  the  production  of  complete  stress- strain  curves 
for  rock  in  compression.  Wawersik  (2)  found  that  visible  cracks  developed 
in  a  marble  in  uniaxial  compression  at  about  90%  of  the  ultimate  load. 

At  a  displacement  about  20%  greater  than  the  peak  of  the  curve,  the 
marble  was  capable  of  resisting  about  half  of  its  peak  load.  Although 
many  large  cracks  appeared,  a  major  portion  of  the  specimen  was  intact 
when  it  was  unloaded  from  this  point  and  removed  from  the  machine. 

Some  rocks,  Wawersik* s  type  II  rocks,  e.  g. ,  are  indeed  incapable  of 
sustaining  load  beyond  their  peak.  These  rocks  have  a  numerically  posi¬ 
tive  "post-failure"  slope  so  that  self-sustaining  failure  can  develop  from 
only  the  elastic  energy  stored  in  the  specimen  and  without  additional 
energy  from  the  machine.  To  test  such  rocks,  it  is  necessary  to  remove 
energy  from  the  rock  in  a  controlled  way  after  the  peak  is  passed.  In 
this  limited  case,  rock  fails  in  a  weakest  link  way;  however,  it  appears 
that  the  behavior  is  very  dependent  on  loading  conditions  (3),  and  except 
as  a  nice  explanation  for  rock  bursts,  it  may  not  be  applicable  to  rock 
structures  at  all.  It  cannot  be  generalized,  then,  that  rock  is  a  weakest 
link  material  in  uniaxial  compression,  it  should  be  obvious  that  in  con¬ 
fined  compression,  rock  sustains  more  load  at  greater  displace n;ent  after 
the  peak  and  is  even  less  a  weakest  link  material. 

What  about  tension?  Hudson,  Brown  and  Rummel  (4)  have  shown 
that  discs  and  rings  loaded  in  diametral  compression  (both  considered  to 
be  indirect  tensile  tests)  also  produce  a  load-displacement  curve  in  which 
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load  is  sustained  at  displacements  beyond  the  peak.  Again,  visible  cracks 
appeared  while  the  specimen  remained  essentially  intact.  Again,  one 
rock  was  found  (Solenhofen  Limestone)  for  which  a  complete  curve  could 
not  be  produced,  but  a  weakest  link  model  was  not  generally  applicable. 

Direct  tension  remains  something  of  a  mystery.  Hughes  and 
Chapman  (5)  showed  a  stable,  bell-shaped  curve  for  concrete,  but  to 
date,  although  there  has  been  at  least  one  serious  attempt  (6),  no  one 
has  succeeded  in  publishing  a  similar  result  for  rock.  On  the  other  hand, 
Brown  and  Singh  (7)  have  shown  that  considerable  acoustic  energy  is  re¬ 
leased  at  loads  considerably  below  ultimate.  Wawersik  (2)  found  loosen¬ 
ing  of  and  fracture  along  grain  boundaries  well  below  the  ultimate  load 
and  concluded  that  "Rock  failure  in  tension  can  be  a  progressive  process 
and  does  not  necessarily  occur  suddenly.  " 

Even  if  we  were  to  concede  that  indirect  tension  of  a  weakest  link 
model  is  applicable,  Hudson  (8)  has  pointed  out  that  direct  tensile  testing 
may  not  be  applicable  to  the  design  of  structures  since  in  practice  grad¬ 
ients  will  always  be  present. 

The  problem  then,  simply  stated,  is  to  make  a  quantitative  model 
for  rock  which  will  allow  progressive  failure  and  residual  strength  to  be 
interpreted  in  laboratory  tests  and  then  utilized  in  structural  design. 

— "Fur,  "  said  the  Deacon,  "t's  might  plain 
Thut  the  weakes'  places  mus'  stan'  the  strain; 

'N*  the  way  t'  fix  it,  uz  I  maintain.  Is  only  Jest 
T*  make  that  place  uz  strong  uz  the  rest.  " 

The  Deacon's  Masterpiece 
Oliver  Wendell  Holmes 

Inhomogenelty  of  Rock 

The  reason  that  the  wonderful  one-hoss  shay  went  to  pieces  "all  at 
once  and  nothing  first,  "  Holmes  tells  us,  was  because  of  the  deacon's 
rather  exacting  construction  standards.  He  recognized  the  weakest  link 


concept  and  eliminated  the  problem  by  making  each  piece  ”uz  strong  uz 
the  rest.  "  One  of  the  most  obvious  things  about  rock  is  that  each  piece 
is  not  "uz  strong  uz  the  rest.  "  This  is  inhomogeneity. 


4 


The  solution  of  most  rock  mechanics  problems  begins  with  a  set  of 
simplifying  assumptions  which  are  nearly  universally  acknowledged  to  be 
invalid.  These  assumptions  usually  include,  but  are  not  necessarily 
limited  to  elasticity  (linear,  small  displacement  elasticity  up  to  failure), 
isotropy  (two  constant  elasticity),  time  independence,  and  homogeneity. 

As  mathematical  tools  have  become  more  generally  available,  these 
restrictive  assumptions  have  been  relaxed  and  we  see  today  a  hugh  volume 
of  literature  treating  rock  as  viscoelastic  (especially  valid  for  salt  and 
potash),  perfectly  plastic  or  work  hardening  plastic  (a  good  model  for 
high  confining  pressure  conditions  as  in  deep  sediments),  transversely 
isotropic  (useful  for  layered  rocks),  etc.,  etc.  Inhomogeneity,  however, 
has  been  largely  Ignored  due  to  .is  severe  insusceptibilicy  to  rigorous 
mathematical  attack. 

Since  the  term  inhomogeneity  seems  to  mean  different  things  to 
different  people,  let  us  clirify  what  will  be  meant  within  the  context  of 
this  study.  By  homogeneous,  we  mean  that  every  point  in  a  material,  or 
at  least  in  a  particular  specimen  of  that  material,  has  the  same  properties 
as  every  other  point  We  must  restrict  the  term  "point"  to  the  continuum 
mechanics  sense  and  ignore  the  atomistic  structure  of  the  material.  This 
condition  of  sameness  of  properties  can  be  violated  in  several  ways. 

The  most  obvious  kind  of  inhomogeneity  is  found  when  the  speci¬ 
men  or  region  of  consideration  consists  of  two  or  more  subregions  with 
distinct  properties.  The  example  that  comes  most  quickly  to  ml  nd  is  a 
sandstone  layer  overlying  a  coal  seam.  A  problem  in  this  situation  is 
deterministic  in  the  same  sense  that  any  homogeneous  rock  mechanics 
problem  is  deterministic;  the  boundary  conditions  and  material  properties 
are  or  can  be  as  well  defined  as  in  the  homogeneous  case.  We  will  call 
this  multiproperty,  multiregion  inhomogeneity  deterministic  inhomogeneity. 
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The  second  way  in  which  the  homogeneous  sameness  of  properties 
condition  can  be  violated  occurs  in  a  rock  type  (often  thought  of  as  a 
partic  'lar  material)  such  as  a  conglomerate  or  large-grained  granite,  which 
is  made  up  of  two  or  more  mineral  constituents.  The  most  obvious  feature 
of  such  a  rock,  even  to  a  layman,  is  its  inhomogeneity.  Unless  we  are 
forced  by  a  particular  test,  such  as  an  indentation  test,  to  sample  a 
very  small  volume  of  the  rock,  its  inhomogeneity  is  usually  shrugged  off 
by  calling  it  macroscopically  homogeneous  and  assigning  single  values 
to  its  elastic  and  strength  properties.  For  many  purposes  this  is  adequate; 
we  will  examine  some  of  the  consequences  later.  We  will  call  this  kind 
of  inhomogeneity  compositional  inhomogeneity. 

A  more  subtle  violation  of  the  sameness  condition  takes  place  in 
rocks  composed  of  a  single  mineral,  such  as  sandstone  or  rock  salt. 

Here,  properties  vary  from  point  to  point,  in  spite  of  mineralogical  same¬ 
ness,  because  of  the  granular  nature  of  the  ;ock.  Each  grain  is  different. 
Each  has  its  own  unique  orientation,  which  may  or  may  not  have  a  pre¬ 
ferred  direction.  Each  may  contain  flaws,  such  as  cleavage  or  fracture, 
or  impurities  in  varying  amounts.  Of  equal  importance  is  the  difference 
in  properties  between  the  mineral  grains  and  their  boundaries.  In  our 
first  course  in  geology,  we  are  taught  to  differentiate  between  quartzite 
end  sandstone  on  the  basis  of  differing  relative  strengths  of  grain  and 
grain  Boundary.  How  then  can  we  consider  those  rocks  homogeneous? 

Good  engineering  results  to  many  problems  can  be  had,  however,  by  tak¬ 
ing  average  or  gross  properties  as  before.  This  kind  of  inhomogeneity 
we  label  granular  inhomogeneity. 

There  is,  of  course,  no  real  separation  between  compositional 
and  granular  inhomogeneity,  with  most  rocks  possessing  some  degree  of 
both.  The  important  point  is  that  this  inhomogeneity  is  not  deterministic. 

In  a  given  problem,  the  boundary  conditions  and  material  properties  are 
not  well  defined.  There  exists  in  such  problems  a  degree  of  randomness. 
We  will  call  this  kind  of  randomness  statistical  inhomogenelty.  K  is 
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with  the  consequences  end  modelling  of  statistical  inhomogeneity  that 
we  will  be  concerned. 

The  discussion  of  compositional  and  granular  inhomogeneity  might 
seem  to  imply  that  statistical  inhomogeneity  exists  only  on  the  scale  of 
the  grain  size.  This  is  useful  in  studying  laboratory  tests  but  becomes 
incomprehensible  at  the  field  scale.  At  this  scale,  deterministic  inhomo- 
gencity  might  include  well  defined  stratigraphy  and  such  structural  features 
os  faults,  dikes,  etc. ,  or  a  particularly  well  mapped  joint  set.  On  the 
other  hand,  a  joint  set  in  which  direction,  spacing,  and  extent  are  known 
only  approximately,  as  is  the  usual  case,  is  an  example  of  statistical 
lnhomogcncity.  Even  within  a  single  bed  or  mass  of  rock,  a  region  may 
contain  more  or  fewer  flaws,  more  or  fewer  inclusions,  or  have  a  different 
mlneralogical  makeup  from  another  region  of  the  same  bed  or  mass.  Such 
variations  may  have  trends  which  can  possibly  be  accounted  for  in  the 
model,  but  this  variation  can,  in  general,  be  said  to  make  rock  statistically 
inhomogeneous  at  the  field  scale. 

Variability  of  Rock  Properties 

Virtually  all  researchers  in  rock  mechanics  are  sufficiently  aware 
of  the  statistical  nature  of  rock  to  avoid  drawing  conclusions  from  the 
results  of  a  single  test.  It  can  be  argued,  however,  that  the  mean  of  a 
"sufficient"  number  of  tests  can  bo  assigned  as  a  "physical  properry"  of 
rock.  Since  a  mass  of  rock  can  be  thought  of  as  made  up  of  an  infinite 
number  of  such  samples,  its  mean  properties  should  be  reasonably  close 
to  the  mean  of  the  sample  properties  and  thus  such  properties  we  used  in 
the  design  of  structures  in  the  rock  mass.  Such  a  practice  obviously  re¬ 
quires  the  use  of  a  sizeable  "safety  factor.  "  The  difficulties  involved  in 
the  use  of  safety  factors  will  receive  extended  consideration  in  the 
section  on  reliability. 

Many  physical  properties,  including  elastic  constants,  unit  weight, 
porosity,  permeability,  thermal  properties,  etc. ,  etc. ,  h  /e  importance 
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In  rock  mechanics  problems;  all  can  be  said  to  vary  to  greater  or  lesser 
extent  from  point  to  point  within  a  rock  mass. 

Hudson  (8)  speaks  of  tensile  suengL.  variation  being  of  three  types: 

(a)  Variation  with  different  specimens.  If  the  same  test  is  performed 
on  the  same  rock  under  the  same  conditions,  different  values  of 
tensile  strength  result. 

(b)  Variation  wit’i  different  specimen  volumes.  If  the  same  test  is 
performed  on  the  same  rock  with  specimens  of  different  volwmes, 
the  average  value  of  the  tensile  strength  changes  with  volume. 

(c)  Variation  with  different  tests.  The  average  tensile  strength  varies 
with  the  typo  of  test  used.  For  example,  the  modulus  of  rupture 
can  be  twice  the  straight  pull  tensile  strength. 

While  Hudson  is  referring  only  to  tensile  strength  tests,  his 
categories  of  variation  are  oqually  applicable  to  otner  strength  measures. 
Wlthin-tcst  variation  is  a  fact  of  all  strength  testing.  We  can  expect 
considerably  different  average  values  for  shoar  strengths  determined  by 
single  shear,  double  shear  and  torsion  tests.  Brown,  Hudson,  Hardy, 
and  Fairhurst  (3)  have  recently  shown  that  between-test  variations  can 
be  expected  with  differing  comprosslve  strength  tests  as  well. 

Let  us  return  now  to  the  problem  of  using  mean  values  from  pro¬ 
perty  tests  as  design  values  for  strength.  Within-test  variations  create 
a  poorly  recognized,  but  very  serious  problem  in  using  test  values  for 
design.  Given  a  weakest  link  model  for  rock  failure,  one  should  clearly 
take  the  mlnimur  of  a  series  of  tost  results  rather  than  the  mean  as  a 
design  value.  Even  without  accepting  weakest  link  failure,  it  would  seem 
prudent  to  design  for  the  worst  condition  likely  to  occur.  This  raiser  the 
question  of  whether  or  not  rock  has  a  minimum  strength  (other  than  zero). 
Particularly  in  tension  it  would  seem  likely  that  if  sufficient  specimens 
were  tested,  one  would  fall  without  application  of  any  load  at  alL  This, 
in  turn,  opens  the  much  larger  question  of  the  sampling  process  Itself. 

In  selecting  blocks  of  rock  from  the  mass,  coring  specimens,  sawing  and 
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grinding  them,  we  reject  at  each  stage  those  samples  which  contain  flaws. 
The  resulting  spoclmcns  can  hardly  be  considered  representative  of  the 
rock  mass.  It  is  for  tho  foregoing  reasons  that  rock  is  often  considered 
to  havo  no  tensile  strength  at  all  in  design  problems. 

Tho  problem  of  a  large  number  of  specimens  giving  a  zero  minimum 
tensilo  strength  is,  of  course,  directly  related  to  the  size  effect.  The 
oxlstoncc  of  a  strength- size  offcct  has  been  widely  debated,  ana  recent 
evidoncc  (9)  seems  to  indicate  tnat  what  appears  to  l>e  <>  si.  e  effect  is 
attributable  instead  to  stress  gradient  effects  or  perhaps  to  differing  avail¬ 
able  onergy  of  tho  testing  machine  in  relation  to  energy  required  for  failure 
of  a  given  volumo.  If  a  sizo  effect  docs  exist,  it  probably  is  a  nonlinear 
docroase  in  strength  with  increasing  specimen  volumo.  if  tho  decrease 
persists  to  largo  volumes,  dosign  values  of  strength  for  rock  structures 
must  be  correspondingly  reduced.  To  apply  laboratory  observed  size 
offocts  to  large  rock  structures  requires  extrapolation  over  several  ordets 
of  magnitude,  wlrlclt  is  always  a  questionable  procedure. 

One  alternative  appoars  to  be  in-situ  strength  testing  involving 
very  largo  specimen  volumes.  Disproportionate  oxj enso  and  the  practical 
difficulties  of  performing  such  tests  have  limited  their  use. 

Another  alternative,  which  the  progressive  failure  model  described 
herein  is  intended  to  explore,  is  to  make  structural  design  problems  as 
deterministic  as  possible  by  geologic  mapping  and  extensive  laboratory 
scalo  testing  and  to  use  tho  results  of  the  laboratory  testing  to  estimate 
distributions  for  tho  statistical  variations  in  the  properties. 

Bctwccn-test  variations  raise  the  question  of  what  test  to  use  to 
obtain  physical  properties.  Particularly  in  tensile  testing  there  is  a  con¬ 
tinuing  proliferation  of  now  tosts  and  little  agreement  oven  on  how  to 
perform  a  singlo  kind  of  test.  It  would  seem  that  less  attention  should 
be  paid  to  the  question  of  which  test  and  how  to  run  it  and  more  to  the 
question  of  why  different  tests  given  different  results.  As  an  example. 
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Hudson.  Brown,  and  Rummcl  0  have  recently  shown  that  the  Brazil 
(diametral  compression)  test,  one  of  the  most  commonly  used  indirect  tests 
for  tensile  strength,  is  not  a  tensile  test  at  all.  With  the  load  surfaces 
in  direct  contact  with  flat  steel  platens  (as  the  test  is  usually  performed), 
they  found  that  failure  always  initiated  in  the  compressive  zone  under¬ 
neath  the  loads.  Even  with  the  load  surfaces  protected,  they  could  not 
cause  failure  to  initiate  at  the  center. 

It  can  be  readily  assumed,  but  perhaps  never  proved,  that  the 
three  types  of  strength  variation  above  arc  manifestations  of  strength 
inhomogeneity.  It  is  much  easier  to  demonstrate  that  variations  in  other 
properties  arc  directly  related  to  statistical  inhomogeneity  of  the  material. 
The  specific  gravity  of  a  rock  specimen,  for  example,  is  clearly  equal 
to  the  weighted  mean  of  tho  specific  gravities  of  all  of  its  subregions 
(if  proper  account  is  taken  of  pore  space).  Similarly,  the  overall  stiffness 
(EA/L)  of  a  uniformly  stressed  specimen  can  be  found  by  a  finite  element 
representation  with  each  clement  given  an  independent  modulus  corres¬ 
ponding  to  the  region  it  represents. 

Probabilistic  Failure  Theory 

All  of  the  above-mentioned  typos  of  strength  variation  can  theore¬ 
tically  be  accounted  for  by  a  model  which  predicts  not  a  unique  failure 
stress  but  a  probability  of  failure  at  any  «trcss  level.  The  best  known  of 
these  models  is  Wcibull's  statistical  theory  of  material  strength  (10,  1 1). 
While  Wcibull's  theory  has  been  shown  to  be  inadequate  for  rock,  probably 
because  it  is  a  weakest  link  theory,  it  will  be  useful  to  examine  it  here. 
The  presentation  here  is  due  to  Hudson  (8)  but  is  available  from  many 
other  sources. 

In  its  simplest  form,  Wcibull's  model  consists  of  V  elements,  each 
having  a  probability  of  failure  Po  at  a  given  stress  o.  The  probability  of 
clement  survival  is  1-Po,  and  tho  probability  of  system  survival  is  (l-PojX 
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or 


If  Pv  is  the  probability  of  system  failure,  then 
1  -  Pv  ,  (1  -  Po)V 

Pv  =  1  -  exp  [V  log  (1  -  Po)]  . 


Wcibull  defines  the  risk  of  failure  R  as  -  V.  (1-Po).  Since  loy 

log 

(1-Po)  depends  only  on  the  stress  j,  the  risk  of  failure;  of  o  differential 
volume  clement  dv  can  be  written  as: 


dR  -  P(s)  dv  . 


If  we  take  R  /  f(c)dv,  then  the  probability  of  the  entire  specimen 
falling  is: 

Pv  i  1  -  cxp[-  /  f(c)  dvj.  (1-1) 

The  problem  here  is  clearly  that  of  defining  Hz)  so  that  the  correct 
cumulative  distribution  of  strengths  (Tv  vs.  failure  load  curves)  will  re¬ 
sult  for  each  test  and  specimen  volume.  The  function  )(  ),  if  properly 
chosen,  would  then  be  a  "true"  material  ptoperty. 

While  the  theory  itself  is  independent  of  the  mechanism  of  failure, 
it  might  bo  thought  of  in  terms  of  a  volume  clement  containing  a  flow 
which  is  "critical"  (in  the  Griffith  sense)  at  the  jmrticular  level  of  normal 
stress  to  which  the  volume  clement  is  exposed.  There  is  no  reason  to 
suppose  that  the  distribution  of  flows  (or  elemental  strengths)  is  well 
described  by  any  commonly  ured  probability  distribution.  Wcibull  assumed 


that  F(o)  was  given  by: 

f(c) 


(1-2) 


whore  c  represents  a  stress  below  which  no  clement  failure  can  occur 
u 

and  c  and  m  arc  material  constants.  He  later  described  it  as  "A  stutis- 
o 

tical  distribution  function  of  wide  applicability,  "  and  indeed  it  is.  This 


function  has  been  used  in  such  diverse  fields  as  econometric  modeling, 
and  the  reliability  of  electron  tubes. 
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Substituting  (2)  into  (1),  we  obtain 

o-o  m 

Pv  1  -  exp  l-  /  [ - J  dv],  (1-3) 

'V  0 

o 

This  is  the  basic  equation  of  Wcibull's  theory.  If  the  threo  con¬ 
stants  arc  truly  material  constants,  and  can  he  found  for  a  given  material, 
Eq.1-3  can  1)0  integrated  to  provide  a  distribution  of  failure  strengths  (in 
terms  of  failure  probability  at  a  given  stress)  for  any  test  or  any  specimen 
volume.  The  process  can  be  quite  involved  since  it  requires  knowing  and 
integrating  the  stress  distribution  ovei  the  entire  spoclmon  volume. 
Nonetheless,  it  has  been  integrated  for  a  number  of  tests.  In  his  first 
paper,  for  example,  Wcibull  presents  the  rosult  that  the  modulus  of  rup¬ 
ture  (3-point  bending  strength)  is  (2m  ♦  2)^m  times  the  direct  pull  tonsilo 
strength.  For  m  =  3,  this  factor  is  2.  0,  and  for  m  -  co,  the  factor  is  1.  0 
(the  strengths  are  the  same). 

Since  its  public  ation  in  1939,  Wcibull’s  theory  has  been  widely 

discussed  and  the  probability  distribution  widely  employed.  It  has, 

however,  been  subjected  to  remarkably  few  experimental  analyses.  In 

many  cases,  experimental  "verification"  has  consisted  simply  of  finding 

a  sot  of  parameters,  c  ,  c  ,  and  m,  to  fit  a  given  set  of  data.  The 

o  u 

existence  of  each  parameter,  of  course,  can  neither  provo  nor  disprove 
the  theory. 

It  would  be  surprising,  in  fact,  if  a  set  of  data  derived  from  any 
statistical  process  could  not  bo  fit  by  a  3-paramoter  model.  The  useful¬ 
ness  of  the  distribution  to  describe  the  fiber  strength  of  Indian  cotton, 
statures  for  adult  males  born  in  tho  British  Isles,  and  the  breadth  of  beans 
of  Phascolus  Vulguris,  lies  in  its  ability  to  "simulate"  a  wide  rango  of 
more  conventional  distributions  such  as  exponential,  normal,  and  log- 
normal,  as  well  as  less  common  shapes. 

Where  the  theory  has  been  tested,  particularly  in  rock,  tho  results 
wore  poor.  Hudson  (8)  found  that  the  theory  would  not  predict  cither  the 
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effect  of  changing  dimensions  or  the  between  test  variations  in  3  and  A 
point  bending  sufficiently  well. 

There  arc  two  entirely  independent  reasons  why  Wei  bull's  theory 
might  not  work.  The  first  lies  in  the  theory  used  in  deriving  equation  1-1. 
Only  if  failure  initiation  in  any  volume  element  tends  to  cataclysmic 
structural  failure  with  no  increase  in  stress  can  equation  1-1  bo  considered 
valid.  V/c  have  previously  seen  that  loev.  is  not  .■  v.i  ■».«  l*:;k  r  .tcn.,1 

under  common  loading  situations.  Equation  I-l  c-mr.ot.  then.  nc  proj>crly 
applied  to  rocks. 

The  second  reason  for  the  failure  of  the  tneory  is  in  the  choice  of 
E(x).  The  function  is  a  material  property  for  any  weakest  link  material. 
There  is  no  justification  at  all  for  selecting  the  V.'eibull  function,  and 
thus  no  reason  to  expect  that  the  constants  are  material  properties.  The 
application  of  equation  1-3  to  any  mate  rial  is  sin. ply  a  curve- fitting 
process,  in  the  Some  sense  as  polynomial  curve  httmn  to  physical  or 
other  data.  It  may  lie  convenient  compel  itionall/  hut  putting  physical 
•inter  pi  elation  oa  the  three  Weibull's  material  constants  :s  as  risky  as 
interpreting  the  coefficients  of  a  curve  fit  polynomial. 

Even  if  wo  reject  Weibull's  theory  as  Inadequate,  it  i^oints  up  at 
least  one*  indispensable*  fact:  in  interpreting  laboratory  strength  test  anJ 
designing  rock  structures,  we  canno'  use  strength  deterministically.  Wo 
must  instead  look  for  and  design  to  a  probability  of  failure. 

Ho  liability 

Quite  naturally,  the  vocabulary  uni  much  of  the  Mathematics  deal¬ 
ing  with  the*  probability  of  failure  arc  loan  1  in  the  literature  of  reli..!  Hr*. y. 
Unfortunately,  workers  in  reliability  arc  usually  concerned  with  the 
probability  of  element  or  system  failure  as  a  (unction  ol  tn  e,  without 
regard  to  whether  it  is  load  or  strength  or  both  which  change  with  time  or 
arc  probabilistic  functions  of  time.  Two  concept;  from  basic  reliability 
theory  are  important  here:  chain  reliability  and  redundancy. 


Chain  reliability  is  represented  by  a  series  circuit  (Fig.  1-1)  in 
which  ouch  element  must  remain  unfailed  for  tho  system  to  be  unfaiicd. 


— -pn - [T| - CSH - GD - E — 

FiG.  1-1  Chain  Reliability 


if  Px  is  tho  probability  that  element  \  is  unfailed,  then  qx  1  -  Px 
is  tho  probability  that  element  x  is  failed.  The  reliability  R  of  tho  system 
is 


R 


p  P  p  p  p 
ABC  DC 


The  probability  of  system  failure  Q  -  1  -  R  is  then: 


Q  • 


1  -  P  P  P  P  P 
A  li  C  D  C 


We  previously  encountered  the  same  notion  in  the  derivation  of 

Weibull's  theory  where  by  taking  the  probability  of  failure  (P  =  1  -  P  -  q) 

as  tho  same  for  all  clemonts,  we  found  the  probability  of  system  failure 
V 

Q  =  1  -  (1  -  P  )  for  V  series  clemonts. 
o 

Redundancy  is  represented  by  a  parallel  circuit  (Fig.  1-2)  giving 
the  opposite  formulation.  Here  it  is  assumed  that  if  any  element  is 
unfaiicd,  the  system  is  unfailed. 


FIG.  1-2  Redundancy 


Taking  qx  os  the  probability  of  element  failure  as  before,  we  have 


and 


R  =  1  -  Q  =  1  -  qAqBqc 


The  resulting  Wei  bull  type  model  for  the  redundant  material  would  have: 

Pv  -  (P  )V  exp  (v  log  (P  )]. 
o  o 

Given  again  that  p  -  qx  is  a  function  of  the  stress,  we  could  write 

Pv  -  exp  (.vfC)  dv).  (1-4* 

We  would  again  b<  faced  with  the  problem  of  finding  f*(  for  the  material 
In  question. 

Notice  that  equation  1-4  describes  a  material  which  loses  r.o 
strength  with  element  failure  while  the  Weibull  weakest  link  material 
lost  all  of  its  strength  with  element  failure.  The  redundant  malctial  is 
clearly  a  |>oorcr  model  for  rock  than  the  weakest  link  material;  the  truth 
is  somewhere  in  bo  tween. 

Anoth  r  appro  ich  to  probability  of  fuller'  ,  whir  h  w  shall  fin*! 
extremely  useful,  is  that  taken  by  K<  cocfoghi  -  nd  (’on  ler  ( 1?).  They 
determine  tire  probability  of  fallen,  of  a  siructuial  con.i^-r.eut  by  as.r.::.  lag 
that  it  is  equal  to  the  probability  that  a  statistically  distributed  siren*, 
exceeds  a  statistically  distributed  strength.  They  enumerate  a  large 
number  of  factors  which  cause  both  stress  and  strength  to  be  probabilistic 
rather  th^n  fixed  quantities.  They  tccommcnd  finding  a  failure  governing 
stress  for  the  component  and  modifying  the  no:  lr.nl  stress  at  the  probable 
failure  point  by  various  factors  using  Monte  Carlo  for  all  variables  con¬ 
cerned  to  produce  an  empirical  frequency  cl  I ;  - 1 .  ibut  ion  for  sties.*..  Sir  ilurly, 
an  empirical  distribution  could  be  constructed  for  i.tmngth  by  starting  with 
the  nominal  (handbook)  strength  of  the  material.  They  say  little  about 
finding  the  distributions  of  the  individual  variables  making  up  the  stress 
and  strength  distributions  but  indicate  that  the  strength  distribution  can 
be  found  by  testing  a  number  of  actual  components  under  the  actual  en¬ 
vironmental  conditions  of  their  use.  This  would  leave  only  the  load 
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distribution  and  certain  environmental  factors  in  the  stress  distribution. 

Kccccioglu  and  Cormier  show  how  to  compute  reliability  R  or  the 
probability  of  comf>oncnt  failure  Q  (  -  1  -  R)  from  the  distributions  of 
strength  (3)  and  stress  (s).  If  f(S)  and  g(s)  are  probability  density  func¬ 
tions  as  shown  in  rig.  J-3,  then  the  probability  that  stress  is  in  the 
interval  ds  about  Sj  is  area  Aj. 

...  ds  ds  .  .  . 

I’fsj  -  —  <  s  <  Sj  +  )  *-  g(sj)  ds  -  Aj 

Tlic  probability  that  S  >  s  is  the  shaded  area  under  the  strength  density 
curve  A.,. 

L  CD 

P(S  >  Sj)  /  f(S)  dS  -  A2 
S1 

The  probability  of  no  failure,  i.  e. ,  the  reliability,  at  is  the  product  of 
these  two  probabilities. 

CD 

dR  -  y(sj)  ds  x  f  f (S)  dS 
S1 

Reliability  of  the  component  is  the  probability  of  the  strength  ex¬ 
ceeding  the  stress  over  the  entire  range  of  applied  stresses: 

CD  CD 

R  =  /  dR  -  /  g(s)  [/  f(S)  dS]ds  (1-5) 

-CD  S 


An  alternative  formulation  can  be  found  from  the  probability  that 
stress  is  less  tl  an  sticngth: 

CD  s 

R  =  /  f(S)  l  /  gfs)  ds]dS  (1-6) 

-  CD  -CD 


and 


By  using  the  properties  of  density  functions: 
CD 

/  f(S)  ds  r  1 

-CD 

CD 

/  g(s)  ds  -  1 

-CD 


3  GENERAL  STRESS  a  STRENGTH 
DISTRIBUTIONS 


it  Cun  presumably  ho  shown  that  1-5  and  I-£  are  equal. 


liquations  i-5  and  1-6  can  be  readily  solved  numerically.  Lxact 
solutions  lor  common  probability  functions  are  obviously  much  more  diffi¬ 
cult.  by  usiruj  an  alternate  approach,  kccecioglu  and  Cormier  give  the 
exact  solution  tor  the  case  where  both  stress  and  strength  are  normally 
distributed 

The  density  functions  for  the  normal  distributions  are: 

y(rJ  ttit  os|j|'  2  1 

s  s 


and 


f(S) 


r  1  ,  s  -  s  . 2 1 

ry  CXpl-T  ( - )  J 

c  27  2  o„ 


1 


where  x  <md  s  are  the  means  and  c  and  o  are  the  standard  deviations. 

u 

Reliability  is  tiie  probability  that  strength  exceeds  stress  or  that 
S  -  b  0.  Taking  S  -  s  -  ",  reliability  is  the  probability  that  >'  >  0.  The 
distribution  of  called  the  difference  distribution,  is  also  normal  and 
Is  thus:  __  2 

hr.)  t-TT-  oxpl  -“(^“H  (1-7) 


C 


where 

and 


S  -  s 

2  ,  2 
°S  +  °s 


Reliability  R  is  given  by  the  positive  density  of  r.  Thus: 

co 

R  :  PC  >  0)  u  /  h(r)dt  (1-8) 

o 

The  probability  of  failure  Q  (  1  -  R)  is  given  by: 

o 

Q  =  P(  <  0)  =  J  h(OdC  (1-9) 

-oo 

This  is  the  shaded  area  in  Figure  1-4.  The  solution  to  equation  1-8 
oi  1-9  can  be  found  from  tables  for  the  standard  noimal  distribution  by  first 
making  the  transformation!:  of  mean  and  standard  deviation  indicated  by 
equation  1-7. 


FIG.L-4  DIFFERENCE  DISTRIBUTION 
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Tne  problem  of  approximating  stress  and  strength  by  normal  distri¬ 
bution**  rap. 'ins.  Without  pie  seating  the  justification  in  detail,  Kecccioglu 
and  CoiinKi  sug.  •  ;t  that  normal  distributions  can  be  approximated  from 
per  ilia  .tie.  i  osl  likely,  and  optimistic  estimates  of  strength  or  stress, 
li  we  take  a.  l».  and  t  as  the  three  estimates  then*. 

a  i  *th  ♦  c 
::  0 

(1-10) 

C  -  a 

‘'s  *  ’  6*~ 

.lit*  considered  good  approximations  They  are  based  on  the  optimistic 
and  pusr.ii  f  lie  e:  tir  ites  being  three  standard  deviations  alx>ve  and  below 
the  n can.  ie  .pectivci/.  The  weights  might  need  to  bo  different  but  the 
concept  is  highly  applicable  to  structural  design  tn  rock  and  Is,  ns  we 
shall  see.  very  much  preferable  to  the  use  of  safety  factors. 

The  safety  far  tor  Is  defined  as  tho  ratio  of  strength  to  stress. 

|ly  taking  t.n<*  ratio  of  tiio  naan  v.ilu  c*f  the  distributed  stress  and  stienglh 
ns  a  safety  factor,  Kecodoglu  and  Cornier  point  out  some  fallacies  in 
the  use  of  safety  factors  for  design. 

Designers  generally  believe  tb.il  o  safety  factor  above  some  pro- 
conceived  value  (for  example,  four  to  eight  for  rock  in  tension)  will  result 
in  no  coin|*c»nont  failure.  On  tho  contrary,  with  these  or  oven  higher 
sifcly  factors,  there  exists  some  finite  probability  of  failure  which  might 
vary  from  acceptably  lov.  to  Intolerably  high.  Similarly,  It  is  commonly 
believed  that  a  safety  factor  of  ono  will  result  in  certain  falluro.  In  fact, 
if  stress  and  strength  distributions  aru  normal  (or  any  symmetric  distribu¬ 
tion).  a  safety  factor  of  one  gives  a  prolcibllity  of  failure  oi  0.  5  regardless 
of  the  standard  deviations  involved.  Even  a  safety  factor  Iocs  than  ono 
does  not  result  in  certainty  of  failure. 

Kccccioglu  and  Cormier  point  out  two  ways  in  which  failure  probab¬ 
ilities  might  vary  with  a  constant  safety  factor. 
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1.  The  mean  of  the  stress  and  strength  distribution  con  bo  changed 
in  the  same  proportion  with  fixed  standard  deviations.  If  both  art*  in¬ 
creased,  the  probability  of  failure  decreases  and  vice  versa. 

2.  The  standard  deviations  can  lx*  altered  while  heaping  the 
means  fixed.  In  this  ease,  an  increase  of  one  o.'  l>oth  standard  deviation*: 
increases  the  probability  of  failure. 

They  quote  examples  where  with  a  constant  safety  factor  o!  2.  i*. 
the  reliability  of  a  component  might  tango  hot  i  .  Of  2d  up  to 

Finally,  it  should  bo  noted  that  Haugen  (Hj  has  written  o  booh 
detailing  ways  for  design  of  mechanical  and  elects ontcchuii! eel  devices  by 
this  method. 


CHAPTER  11 


A  ONE- DIMENSIONAL  MODEL 


I IndSOU's  r/cvlol 

The  chief  criticism  of  the  Woibull  theory,  insofar  us  its  npplice- 
tlon  to  roe!.  is  concerned,  Ik  that  It  Is  a  weakest  link  model.  An  attempt 
to  foimulale  a  Welbull -type  treatment  based  on  a  redundant  model  gave 
an  unacceptable  result  because  it  allowed  for  no  deterioration  with  clement 
failure,  creating  a  sort  of  strongest  link  model. 

A  simple,  mathematically  elegant  one-dimensional  model  which 
allows  for  structural  deterioration  with  element  failure  has  been  con¬ 
structed  by  Hudson  (8).  This  model  nakos  use  o(  distributions  of  slrons 
and  slrongth  In  a  manner  strikingly  similar  to  Kccecloglu  and  Cormier's 
reliability  computation.  There  is  an  important  difference,  however. 

'Rather  tiian  Irclng  concerned  with  probability  of  failure  of  an  entiro  com¬ 
ponent,  Hudson's  model  looks  at  the  piobability  of  failure  of  an  Infinite¬ 
simal  element,  and  deteriorates  the  structure  in  a  systematic  way  as  each 
element  (alls. 

The  model  13  most  easily  visualized  os  a  largo  number  of  parallel 
elements,  similar,  as  Hudson  points  out,  to  the  surings  In  a  mattress. 

Hudson  works  entirely  in  torins  of  strain  rather  than  stress,  and 
further  assumes  that  the  assemblage  is  to  lie  tested  in  a  "stiff"  (strain- 
controlled)  machine. 

As  the  assemblage  is  loaded,  a  mean  strain  (disturbance)  is  applied 
to  the  specimen.  The  strains  on  each  elemont,  however,  arc  assumod  to 
bo  distributed  about  this  mean.  Hudson,  as  did  Kecccioglu  and  Cormier, 
assumes  a  normal  distribution  for  illustration;  tho  model  is  not,  however. 


dependent  on  the  choice  of  distribution,  lie  assumes  that  each  element 
is  linearly  elastic  until  its  strain  reaches  a  certain  value*,  at  which  failuie 
occurs  and  the  element  can  no  longer  sustain  lend.  As  the  specimen  is 
loaded  liy  the*  strain-controlled  machine,  th''  mean  strain  is  presumed  to 
increase  monotonically.  1.  e. ,  it  is  the  independent  variable  and  the  lotnl 
load  or  stress  on  the  assemblage  is  the*  response. 

Hudson  unaccountably  lakes  the  variances  of  me  two  distributions 


**d 


2  ?. 
A  and  r 


.  2  ? 
r  a 


Tito  effect  of  this  assumption  is  that  as  is  increased  through  the  test, 

the  variance  of  the  *  .  distribution  also  increases,  Hole  that  «•  remains 

d  s 

constant  so  ti  nt  the  assumption  that  variance  is  a  function  of  mean  value 

2  2 

has  no  effect.  Tire  quantity  %T A  is  simply  an  arbitrary  constant. 

s 

Hudson  (Kiints  out  that  x  and  1  should  be  chosen  so  that  Z ,  and  V 

d  i. 

arc  always  greater  than  three  standard  deviations  away  from  zero.  This 
*wJ  11  reduce  the  number  of  negative  strains  and  strengths  to  about  0.  )3,. 
of  the  total.  Unlroundod  distributions  hive  the  slight  disadvant  »ge  of 
always  producing  some  negative  values.  Because  standard  deviation  is 
the  square  root  of  the  variance,  the  requirement  is  satisfied  by  x  end  1 
themselves  exceeding  three. 

As  is  increased,  elements  fail  in  succession;  failure  occurs 
in  each  clomont  when  <*.  At  any  value  of  V^,  the  probability  of  fail¬ 
ure  is  given  by  the  difference  norma!  distribution  as  shown  by  KccccJoglu 
and  Cormier. 

Hudson  next  makes  the  crucl<il  assumption  that  the  resist  nice  of 
the  overall  assemblage  Is  proportional  to  I  (C^) ,  the  "survival  function" 
or  i*ui*ortion  of  elements  unfailed  at  any  c Thus; 

«r=  ci  £d> id  oi- 

where  c  -  mean  stress  response 

and  C  a  constant  of  proportionality. 


As  «*  approaches  zero,  the  survival  fraction  I  approaches  one, 
thus  the  constant  C  is  the  tangent  modulus  at  zero  load  or  Young's  Modulus 
of  the  material.  The  quantity  C  1  is  the  relation  between  stress  and  strain 
at  any  r*c*lnl  and  is  therefore  the  secant  modulus. 

Resulting  theoretical  stress-strain  curvos  for  the  assumption  that 
h  1  arc  shown  in  figure  11-1.  In  general  shape,  these  curves  bear  a 
marked  similarity  to  complete  stress -strain  curves  for  rock  In  compression. 
Koto  also  that  for  I:  1  -  co,  tire  "perfectly  homogeneous"  case,  the  stress 
drops  lie  lantarK'ou..ly  to  zero  and  the  model  represents  a  weakest  link 
material. 

*iho  behavior  shown  in  figure  11-1  can  bo:.t  bo  described  as  "pro¬ 
gressive  failure*  or,  as  Hudson  pats  it,  "structural  breakdown.  "  Clearly, 
such  a  procoss  represents  an  apt  way  to  describe  the  way  rock  falls. 

Development  of  a  Numerical  Model 

In  older  to  investigate  some  of  the  ramifications  of  Hudson's  modol, 
Y»  discrete  numerical  model  v/us  constructed  and  prog r*»tr. mod  for  the  digital 
computer. 

There  ore  o  finite  number  of  elements  N  arranged  in  parallel  as 
before.  It  is  customary  in  rock  mecltenics  work  to  use  original  specimen 
dimensions  in  computing  stross  and  strain.  The  increase  and  later  reduc¬ 
tion  of  oron  in  a  compression  specimon  is  ignored  and  strain  is  tuken  as 
‘•engineering*  rather  than  "true"  strain.  The  modol  can  bo  simplified, 
therefore,  v/lthout  sacrificing  generality,  by  working  in  tcrni3  of  force 
and  displacement  rather  than  stress  and  strain. 

Hudson  assumes  that  the  reduction  in  stiffness  as  each  element 
fails  is  a  constant.  The  total  silliness  K  of  a  collection  of  parallel  elastic 
elements  is  the  sutn  oi  element  stiffnosscs  Y.  ,  thus: 

N 

K  -  >  k  NK  . 

1  1  1 
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Tii*  lolril  force  I’  on  the  assemblage  Is: 

N  N 

r  =  :  f  -  :  k  r  r  n  r  k  01-2) 

i-i  1-1 

where  -  the  distributed  element  displacement 

end  t'f  the  mean  element  displacement. 

As  the  distribution  of  Is  incensed,  p.ofMgaling  the  distuib.  nco 
into  the  assemblage,  the  individua  l  will  begin  tc*  exceed  the  distributed 
strengths  r.  .  We  will  assume  th.it  they  lose  t'p  »r  lojd  hociiipy  t  ikicHv 
ns  before  and  that  we  should  deal  only  with  the  number  of  uufoiled  ele¬ 
ments  n  whuro: 

n  =  1  N  (11-3) 

and  I  Is  tho  survival  fraction. 

Combining  equations  11-2  and  11-3,  we  get: 

I  -  I  ni.  7  ■-  I  K  T  (11-4) 

tvhich  is  Hudson's  equation  (II- 1)  in  tonus  of  foicc  end  displacement. 

Koto  that  I  is  again  a  function  of  the  difference  distribution,  this  time 
of  (C(  -  c4).  In  equation  11-4,  K  is  tho  initial  stiffness  or  tangent  stiff¬ 
ness  and  1  K  Is  tho  secant  stiffness. 

l/’t  us  first  ask  what  would  happen  If  wo  were  to  distribute  forces 
rather  than  displacements.  The  element  displacements  would  nit  11  bo  dis¬ 
tributed  In  the  sumo  way  since  f^  -  i:  and  equation  (JI-4)  would  still 
ensue.  The  survival  fraction  would  still  be  the  result  of  tho  difference 
distribution  and  the  forr  of  the  resulting  stress-strain  curves  would  be 
unaltered. 

The  result  in  rather  surprising  since  tho  form  of  tho  complcto  stress- 
strain  curve  is  lather  different  in  a  laid -con trolled  (soft)  machlno.  Wo 
see,  however,  that  working  In  terms  of  forco  has  not  modeled  soft  loading 
at  all  since  it  Is  1^  that  If.  increased  monotonically  and  not  r. 
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To  make  F  the  independent  variable,  v/c  need  a  rule  for  distributing 

n 

l.  among  the  surviving  elements  so  that  V  -  ;  f.  .  With  our  previous 
1  i-l  ‘ 


assumption  that  all  k^  wcio  the  same,  the  simplest  rule  is  to  make  all  = 
and  thus  F  =  /.  This  simply  s«ys  that  the*  progressive  failure  of  the  us.tlei- 
ial  results  entirely  from  inhomogeneous  strength. 

Two  other  ways  are  appio-nt.  Tec  first  is  to  keep  J:  k  <nd  to 

assume  that  dltfering  f.'s  result  from  a  ei  tubaLion  of  ir.iti  !  *.  1  cements. 

j 

V/e  must  have  t  0  so  tlu.t  •  d  at  i‘  (>.  Then,  f,  i:(  .  -t  ).  This 
o  i  i  o 

could  also  be  icpiescntcd  as  an  iriti  d  force  dicti i’>  ‘.ic.a  f  1:  ■  .  The 

i  o 

assumption};  that  k.  k  and  ol  initial  foices  or  dir.pi  cc*t?.f  nt*;  an*  physically 
less  appealing  than  the  second  possibility.  Hero  we  assume  that  the 
disliihution  of  f  results  from  n  distributed  k  and  that  <>) )  -  /. 

1  J 

The  second  approach  was  chosen  for  proernmninn  to  test  the  clfect 
of  the  distribution  shape  and  r.i;*e  on  the  resulting  sti ess-strain  {actually 
force- displacement)  curves.  Hudson  had  compared  ihc  losulls  of  using 
l»oth  distribution:,  normal  air?  noth  ur.  forr  ird  foun?  that  the  < 'lives  ware 
almost  indistinguishable.  1  hi n  v.ms  thought  lobe*  inconcler.iv*,  lov..  vo\ 
and  the*  five  1  rounded  distributions  discussed  in  A:*j*endi>:  A  were  piogiainmod 
so  that  the  stiffness  and  strength  distributions  could  hr  selected  with  any 
shape,  any  moon,  and  any  width  or  variance,  desired. 

The  model  was  first  programmed  with  a  foice  failure  criterion  and 
arranged  so  that  U  i  a  given  ru  t  of  (joncrated  by  Monte  Cailo  horn 

the  selected  distributions,  both  iorcc  an  1  displacement  loading  icsnlts 
wore  idotted.  Since  tlic  compatibility  equation  cays  that  all  ^  /,  the 

force  fj  on  an  element  is  ^  1;^  holme  occ  us  when  «•,  ,  the 

strength  oi  the  element  (now  in  force  tcims).  Ti»o  o-dor  of  failure  is  thus 
determined  hy  the  ratio  s^/k,  which  is  actually  the  failure  value  of  f  . 

We  sec  that  the  supposed  force  f  lihnc  niterion  has  Income  a  displacement 


failure  criterion  duo  to  the  restrictions  of  the  model.  This  approach  was 
tried  with  Interesting  results,  but  it  was  discovered  that  the  rotio  of  two 
random  variates,  even  from  tho  simple  uniform  distribution,  has  a  complex 


M  c 
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distribution  which  indites  pliysic.il  interpretation  of  the  effects  of  distri¬ 
bution  shape  i ikI  si/c  i.ithei  difficult.  A  second  vorsicn  of  the  program 
was  written  in  which  the  displacement  failure  criterion  was  generated 
directly. 

Let  us  now  look  at  the  workings  of  this  program.  The  number  of 
elements  is  lead  in  along  with  the  shape  and  variance  of  the  distributions. 
Subroutines  arc  called  which  gen. rate  a  stiffness  and  a  s!  length  value 
(now  in  dl.sp!  icon  cut  units)  for  each  element.  The  elements  are  then 
soiled  and  rerun  bored  in  order  of  increasing  fail uro  displacement  and 
the  total  stiffness  of  the  assemblage  Is  computed.  In  displacement  load¬ 
ing,  the  total  force  is  calculated  from  the  given  displacement  and  the 
total  stiffness;  the  stiffness  of  the  failed  element  is  subtracted  to  obtain 
the  new  total  stiffness  and  the  pioccss  is  repeated.  Tin*  arrays  contain¬ 
ing  total  force  and  corresponding  displacements  are  saved  and  can  bo 
printed  out  or  plotted  or  both. 

loco  loading  requires  only  t  sltghl  modi  beat  ion  of  tiro  displace¬ 
ment  procedure.  At  each  step,  a  check  must  bo  made  to  see  if  the  failure 
of  one  element  has  redistributed  forces  in  such  o  way  as  to  cause  failure 
of  additional  elements.  With  a  force  failure  criterion  this  can  be  done 
by  dll  eel  comparison  of  new  element  force  to  strength.  It  is  much  simpler 
(and  more  easily  adapted  to  displacement  failure),  however,  to  check 
for  a  drop  in  total  force  at  each  step.  If  such  a  drop  occurs  in  the  dis¬ 
placement  loading,  it  is  ignored  in  force  loading  and  the  previous  force 
substituted. 

Since  drops  in  force  are  always  small  and  localised  phenomena  on 
the  rising  side  of  the  force  displacement  curve.,  the  two  methods  tend  to 
give  nearly  identical  plots  up  ;o  the  peak  of  the  curve.  After  the  peak  Js 
passed,  the  force  loading  curve  merely  displaces  without  limit  at  a  con¬ 
stant  force  and  is  uninteresting.  After  some  testing  of  these  fdets,  the 
computing  and  plotting  of  the  force  loaded  curve  were  discontinued.  The 
final  version  of  the  one-dimensional  program  is  presented  in  Appendix  B. 
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We  note  that  the  plotted  path  does  not  agree  in  detail  with  either 
the  force  loading  path  or  the  displacement  loading  path.  The  drops  in 
force  in  the  displacement  loading  path  become  very  small  unless  the  drop 
in  stiffness  is  a  large  pioportion  of  the  remaining  total  stiffness  os  is 
grossly  exaggerated  in  these  examples.  Deviations  (if  1  ho  plotted  path 
from  the  displacement  loading  path  are  I  hero  fore  negligible  in  large 
assemblages  of  elements.  Up  to  the  pea':  of  the  cure  o: . J >  mil  r  yener.J 
drops  in  force  will  cause  large  deviation.;  o'  the  fore  »  'coding  path  from 
the  plot.  l.'xecpt  where  such  a  general  drop  might  occ‘'r,  it  can  be  claimed 
that  Ihj  plotted  output  fairly  represents  either  force  01  displacement  load¬ 
ing  up  to  the  peak  of  the  curve. 

The  model,  as  developed  above,  can  theoretically  represent  any 
rock  loading  situation  In  which  stress  gradient:;  can  be  assumed  negligible. 
This  essentially  precludes  modeling  all  but  uniaxial  tension  and  com¬ 
pression  tests. 

To  reduce  the  number  of  independent  variables  in  the  model  and 
facilitate  comparison  between  tests,  a  hypothetical  rock  was  selected  fo» 
modeling  in  a  uniaxial  compression  test.  The  specimen  is  a  cylinder 
with  a  diameter  of  2  inches  and  a  length  of  •!  inches.  The  rock  has  n 
modulus  of  7.  6  x  10**  psi  giving  a  stiffness  K  ~  of  (>  x  10^  lbs.  It  lias 
n  yield  point  of  7600  psi  or  24,  000  lbs,  and  It  has  a  strength  of  about 
13,  000  psi  or  40,  000  lbs. 

Results 

The  ono-dimcnsional  model,  as  programmed  in  Appendix.  R,  has  a 
large  number  of  independent  variables.  Clearly,  a  factorial  design  in¬ 
volving  sufficient  levels  of  all  variables  is  out  of  the  question;  a  system¬ 
atic  reduction  was  attempted  instead. 

The  logical  starting  point  Is  the  number  of  elements  N.  It  seems 
reasonable  to  assume  that  the  effects  of  varying  N  will  appear  independently 
of  the  other  variables  in  the  model  (i.  o. ,  t ho  distributions).  These  effects 


dio  of  two  distinct  kinds.  The  least  Important,  but  probably  most  obvious 
in  the  f i c.ui v'ii  which  follow,  is  that  the  plotting  of  smaller  numbers  of 
points  by  a  "connect  the  dot.;"  rulo  will  produce  a  jagged  curve.  More 
important  j  -  the  statistical  "law  of  laigc*  numbers"  which  says  that  if  a 
"sufficient"  number  of  samples  Is  taken,  the  sample  distribution  is  arbi¬ 
trarily  <  lo.;o  to  the  theoretical  distribution.  This  moans  that  as  N 
approaches;  infinity,  the  curve  produced  by  the  model  bocomes  determinate. 
The  ntmibet  of  el  aments  should  bo  selected  large  enough  to  make  the  farce 
displacement  curve  smooth  and  reproducible. 

This  requirement  must  be  balanced  against  the  increased  cost,  in 
torms  of  computing  *imo,  of  increasing  N.  Figures  11-2  through  11-6  show 
the  results  of  vatying  N  (N  -  10,  100,  500, 1000)  with  both  distributions 
uniform.  The  curves  for  N  -  500  and  N  -  1000  aro  smooth  and  noarly  in¬ 
distinguishable,  whllo  the  curves  for  N  =  10  and  N  -  100  are  rough  and 
not  reproducible.  It  was  concludod  that  500  elements  were  sufficient, 
and  this  number  was  used  in  meet  of  the  later  trials. 

Originally,  testing  was  done  with  A  ,  the  lower  bound  of  the 

s 

strength  distribution,  fixed  at .  004  inch.  This  was  to  simulate  a  rook 
with  a  yield  point  of  24,  000  lbs.  or  about  7, 600  pnl.  The  modol  Is  elastic 
until  the  first  element  fails  which  will  occur  at  the  displacement  where 
tho  I  oast  strength  Is  exceeded.  The  stiffness  returns  to  zero  when  the 
last  elomont  falls,  thus  B^  represents  a  "maximum"  displacement  The 
result  of  changing  the  variance  of  tho  strength  distribution  is  shown  in 
Figure  11-7.  Hero  the  stiffness  distribution  was  held  constant  while  the 
variance  of  the  strength  distribution  ranged  from  xoro  to  5. 6  x  10 
Both  distributions  wore  uniform  for  this  tost  With  bounded  distributions 
the  role  of  the  strength  Vdrlanco,  or  specifically  tho  upper  and  lower 
bounds  of  tho  strength  distribution,  is  to  fix  the  yield  point  and  tho  max¬ 
imum  displacement  at  which  the  specimen  can  bear  load.  If  A  and  B 
arc  coincident  (variance  -  0),  tho  yield,  peak  and  maximum  displacement 
occur  at  the  same  displacement  and  weakest  link  failure  occurs. 
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The  effect  of  V,,  it  shown  similarly  in  l'lg.  11-8.  Here  V  ranges 
from  JMsro  to  10  with  both  distributions  .gain  uniform.  Th*  curves  are 
quite  intcrclwingoable.  This  unexpected  result  is  due:  to  the  Jndcpea 
of  the  sllflm  js  oik!  strength  dishlbutionr..  In  ihe  progt  c  tiu.  site  r.jths 
are  sorted  in  order,  but  because  it  is  inJo{*cn:!cnt  l!io  stil'iu  n  .  is  Milt 
quite  random.  Over  u  largo  number  of  .imaf!  «JispI».cei.-oM  .‘-.tops,  ther 
foro,  the  drop  in  stlUnc:>&  at  each  i  lep  will  lend  to  the  i  can  of  the  « •  i r. - 
UibeUon  which  is  determined  bcJcithjiui  by  select.?',  j  io;  a  to*al 
stiffnes.  of  i'  x  10^  lbs,  each  element  l»as  n  mean  stiff  ”  of  t  10^/K 
lbs.  The  fact  that  the  force  displacon  ;nt  curve  is  Indc'.vndvnt  of  the 
•tiffnoss  distribution  is  again  Illustrated  in  Fig.  11-9  where  thr.-  strength 
distribution  is  again  uniform  and  fixed  and  the  stiflness  distribution  shape 
it  soli  is  varied.  Again  tiro  curves  are  interchangeable. 

Figures  11-10  through  11-14  show  the  effect  of  the  flvo  jjonsiblo 
shapes  of  the  strength  distribution  on  the  force— displacement  curve. 

Figure  11-10  shows  the  smooth,  rathe?  sy:r.i.;ctr:cu)  curve  produced  ny  the 
uniform  distribution.  Fig.  11-11,  with  tho  strength  unimode lly  ciisli  Touted, 
is  still  smooth  but  has  a  peak  about  one-third  higher  than  Fig.  11-10. 

This  is  because  fewer  elements  fail  In  tho  oarly  stages  of  displacement 
and  tho  sporfmon  retains  more  stiffness  et  higher  displacements,  rig  11-12 
•hows  a  still  higher  poa):  caused  by  oven  less  element  failure  at  small 
displacements  with  tho  left  skewed  distribution.  Fig.  11-13  shows  tho 
opposite  effect  of  the  right  skewed  distribution  where  elements  (all  most 
rapidly  In  cnrly  loading  and  slower  near  the  end,  thus  flattening  the  full- 
iri'i  side  of  the  face  displacement  curve.  Finally,  Fig.  11-14  shows  the 
behavior  when  tho  drstrii>ution  is  bimodal,  causing  a  plateau  effect  in 
tho  middle  range  where?  few  elements  arc  falling. 

Up  to  this  point,  the  minimum  of  tho  strength  distribution  A  was 

s 

fixed  at  0.  00*1  inch.  In  Fig.  11-15, Its  effect  is  analyzed  Independently  with 

the  strength  distribution  uniform  and  the  minimum  strength  A  varying  bom 

R 

0  to  A.  008  inch  and  the  strength  variance  huh1  constant. 
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FIGURE  11-14  Strength  Distribution  Bircodal 
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>»*.•;  mo del  used  a  fixed  mean  rathci  than  a  fixed  minimum 
tin  tit.  a  it  ;lii  d i . j •  t  * '  ITi;-;.  Il-iu  through  11-20  show  that  the 

mirneu;  .1  i"u.i  *1  }/.  t.H-e.;  the  characteristic  "Hudson*'  family  of  curve;.  Jf 
the  nn-  at  i.\  fixed,  in  this  case  to  U,  Ul5  inch. 

liegurdles.-  t>‘  V .  varying  variance,  each  family  of  curves  is  soon 
to  rdtairi  common  |.vinL.  Tin's  j  oint  is  always  located  at  the  displacement 
corresponding  tu  th  •  u  :an  of  tit.  stieiigth  distfibution.  i'or  the  three 
symmetric,:]  c’isln’bu'.  ions  (unite  rr.,  uiumod.il  and  bi modal),  ti  e  force  at 
th*.  sh  .o\!  point  In.  the  pe.ih  force  for  the  homogeneous  (no  variance)  cur vo. 
This  is  because  at  the  mean  sU»  nglh  half  of  the  elements  have  failed  nr.d 
the  stiHn  is  thus  cut  in  hall.  Similarly,  the  skewed  curves  have 
shared  points  at  one-third  and  two-thirds  the  peak  force  of  the  homogen¬ 
eous  curves. 

l'inally,  we  see  that  the  curves  for  the  no  variance  case  arc  all 
identical.  This  unrealistic  behavior  is  tho  only  possible  result  of  assum¬ 
ing  that  rock  is  uniform  and  tha!  stresses  are  uniform  in  a  uniaxial  lest. 

Coned  ns  joins. 

Thu  one-dimensional  numetical  model  agrees  qualitatively  with 
Hudson's  results  if  similar  assumptions  are  made.  The  use  of  bounded 
distributions  and  numerical  techniques  allows  the  effect  of  all  variables 
to  be  examined  in  detail.  In  addition  it  allows  physical  interpretations 
to  be  placed  on  parameters  as  follows: 

1)  The  C  in  Hudson's  model  (Eq.  11-1)  is  the  Young's  modulus 
of  the  material. 

2)  C  I  is  tho  secant  modulus  at  any  displacement. 

3)  The  distribution  of  stiffness,  if  independent  of  strength, 
lias  no  effect. 

A)  The  minimum  of  the  strength  distribution  determines  tho 
yield  point. 

3)  The  maximum  of  the  strength  distribution  determines  whoio 
the  force  displucmoni  curve  n  turns  to  the.  axis. 
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FIGURE  11-16  Effect  of  varying  with  constant 

mean  strength 
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(•)  Tin-  sh.ijv  o t  ti.*;  fuuv  displacement  curve  can  be  quali- 
laiivi  lv  <!•  ,;i  imcci  in  terms  of  tin*  shape  of  the  strength 
<Ji r.tr  i’.  utit  a. 

/  it’. oa  -.i  »t  ».m  i  .  t.nf.  .  -:i;  t0  qU  .utttetivoly  model  a  particular 
roe*,  in  vi  uniaxial  tv  ..t  u.tinii  taif*.  method,  t!io  model  h<is  great  value  in 
intoipiet  i, :  tin  pro  ;i osufv^  failure  of  rock  in  terms  of  inhomogenoity  of 
the  in  1 1 v  1 1  i). 
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APPENDIX  A 
DISTRIBUTIONS 

The  five  distribution  functions  used  in  this  study  ere  bared  on  a 
set  of  functions  suggested  by  Hemnerle  04)  for  digital  computing.  The 
first,  the  uniform  distribution,  is  well  known,  and  the  others  arc  derived 
from  it  in  the  following  way: 

1)  A  variate  from  the  uni  mod  a  l  distribution  is  obtained  by  adding 
two  random  variates  from  the  uniform  distribution. 

2)  A  variate  from  either  of  the  skewed  distributions  is  obtained  by 
doubling  either  the  largest  (RSICEW*)  or  the  smallest  (L3KEW*)  of 
two  variates  from  the  uniform  distribution. 

3)  A  variate  from  the  bimodal  distribution  is  obtained  taking  a  variate 
from  the  unimodal  distribution  and  subtracting  from  1  if  less  than 

1  and  from  3  If  greater  than  1. 

Most  large  computers  have  available  a  subroutine  which  produces 
pseudo-random  uniformly  distributed  numbers.  At  UWCC  (the  University 
of  Wisconsin  Computing  Center),  a  function  RAN  UN  which  can  be  called 
from  Fortran  programs,  returns  uniform  variates  in  the  range  0  to  1  (1 5). 
Since  tho  uniform  variates  are  in  the  range  0  to  1,  the  procedures  described 
above  for  the  other  distributions  produce  numbers  in  the  range  0  to  2. 

In  order  to  provide  numbers  in  the  range  a  to  b,  a  scaling  technique  sug¬ 
gested  by  Naylor,  ct  al.  (16),  was  used. 

To  test  the  generation  and  scaling  routines  cumulative  distribution 
functions  were  found  for  1,  000  variates  in  the  range  2  to  8  for  each  dis¬ 
tribution.  These  are  shown  in  Pigs.  A-l  to  A-5  and  are  seen  to  match  very 
well  the  functions  expected  from  integration  of  the  continuous  distribution 
functions. 

Subroutine  RSKEW  actually  produces  a  distribution  which  is  loft  skewed, 

1.  o. ,  the  mean  of  the  distribution  is  to  the  right  of  the  midpoint  with  the 
toil  of  tho  distribution  stretched  out  to  the  left.  Similarly,  LSKLW  produces 
a  right  skewed  distribution. 


figure  (A-I)  CDF  FOR  UNIFRM 
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FIGURE  (A-2)  CDF  FBR 
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figure (a-3) CDF  FOR  LSKEW 
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It  should  be  noted  that  except  for  the  uniform  distribution  the 
functions  used  arc  not  stochastically  derivable  (iuslitinblc)  in  tire  sense 
that  the  normal,  gamma,  and  other  commonly  used  in  simulation  arc. 

There  is  no  ioason  whatever  to  expect  that  these  distributions  represent 
the  properties  of  real  rock.  (There  is,  however,  nothing  which,  at  the 
present  time  at  least,  would  lead  to  on  expectation  that  any  other  distri¬ 
bution  represents  rock  properties.  Tins  problem  is  unstudied  and  relatively 
unstudiahle  as  well).  They  were  chosen  because  they  arc  cheap  to  gen¬ 
erate*  on  n  digital  computer  (a  single  normal  variate  requires  the  gonerai ion 
of  12  uniform  variates  by  the  Teichroew  approximation  (16)  or  horn  6  to 
12  times  longer  than  methods  used  here),  and  because  they  jo  present  a 
sufficient  variety  of  behavior  that  the  effect  of  distribution  shapes  on 
progressive  failure  could  be  thoroughly  examined.  An  added  benefit  is 
that  the  distributions  are  all  bounded  on  a  range  of  a  to  b,  allowing  physi¬ 
cal  interpretation  to  be  placed  on  the  bounds  a  and  b. 

Ac  proposed  by  liommorie,  the  statistics  of  the  \  derived  distribu¬ 
tions  had  not  boon  studied,  and  calculation  of  the  mean  tmd  variance 
wore  required  according  to: 

b 

E(x)  =/  x f(x)  dx  (A-l; 

a 

and  , 

2  k  2 

E(x  /  r  /  x*  f(x)  dx  (A- 2) 

a 

The  results  of  theso  Integrations  are  shown  in  Table  A-l.  Also 
shown  is  the  width  of  each  distribution  (b-a)  in  terms  of  lh°  variance 
V(X)  where  V(X)  io  given  by: 

(A-  3) 


v<x)  :  i:(x2)  -  [i;  (v)]2 


TABLE  A-l 

STATIST; CS  TOR  FIVE  DISTRIBUTIONS 


No. 

Distribution 

E(X) 

E(X?) 

b-a 

1 

l)  niform 
(UNI1RM) 

a  +  b 

o?  4-  ab  4  b^ 

12V 

2 

3 

2 

Unimodnl 

a  +  b 

7a^  4-  10ab  4  7b^ 

12V 

(UNI  MOD) 

2 

24 

3 

Right  Skewed 
(LSKEW) 

b  4-  2a 

3 

b^  4-  2ab  4-  3a^ 

6 

18V 

4 

Left  Skewed 

2b  4  a 

3b^  4  2ab  4  a^ 

18V 

3 

6 

5 

Bimodal 

a  +  b 

2  2 

3a  4  2ab  4  3b 

8V 

2 

8 
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APPENDIX  B 

ONE-DIMENSIONAL  PROGRAM 

The  following  pages  contain  a  listing  of  the  main  program  and 
seven  subroutines,  all  in  Fortran  IV  which  were  used  to  produce  the  force- 
displacement  plots  in  Chapter  II.  The  programs  contain  references  to 
numerous  subroutines  provided  by  the  University  of  Wisconsin  Computing 
Center  (UWCC)  as  system  library  routines.  Most  of  these  arc  involved 
with  the  graphics  package  which  produced  the  plots.  The  operation  of 
the  programs  is  described  briefly  below.  While  the  programs  as  presented 
can  be  run  only  on  the  Univac  1100  at  UWCC,  it  is  hoped  that  some 
insights  into  the  model  can  be  gained  from  their  presentation  here.  It  is 
further  hoped  that  the  distribution  subroutines  will  be  useful  to  others 
interested  in  simulation. 

The  main  program  first  fills  arrays  containing  element  stiffness 
(STirr)  and  strength  (DISM)  with  random  variables  from  the  appropriate 
distribution  by  calls  to  GENS  and  GENIC.  It  computes  initial  total  stiffness 
(TSTD  and  uses  a  UWCC  utility  routine  (URSORT)  to  sort  both  arrays  in 
ascending  order  of  strength.  In  displacement  loading,  the  force  (FORCE) 
at  oach  load  increment  is  computed  as  the  product  of  current  total  stiff¬ 
ness  and  the  failure  displacement  (D1SPL).  The  total  stiffness  is  then 
reduced  by  the  stiffness  of  the  failed  element.  The  arrays  (FORCE  and 
DISP1J  are  then  scaled  and  plotted. 

Subroutines  GENK  and  GENS  fill  arrays  STIFF  and  DISM,  respect¬ 
ively,  with  random  variables  by  calls  to  appropriate  distribution  sub¬ 
routines.  The  so -array- filling  routines  also  compute  sample  statistics 
for  comparison  with  distribution  statistics. 

Tho  distribution  subroutines  (UNIFRM,  UNI  MOD,  RSKF.W,  LSKKW, 
BIMOD)  obtain  random  variates  from  tho  five  distributions  as  described 
in  Appendix  A. 
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MECHANICAL  BEHAVIOR  OF  ROCK  UNDER  CYCLIC  FATIGUE 

1  2 
By  Bezalel  C.  Halm  ton  and  Chin  Man  Kim 


INTRODUCTION 

Rock  and  rock  structures  such  as  open-pit  benches,  mine  esca- 
vations,  bridge  abutments,  dam  and  road  foundations  undergo  cyclic 
loading  caused  by  earthquakes,  traffic,  drilling,  blasting,  etc.  This 
type  of  loading  often  causes  a  material  to  fail  at  a  stress  lower  than 
its  determined  strength,  a  phenomenon  called  fatigue.  Fatigue 
characteristics  are  usually  presented  in  the  form  of  an  S-N  curve 
(stress  level  versus  the  number  of  cycles  required  to  bring  about 
failure).  The  top  stress  value  that  appears  to  be  unaffected  by  cyclic 
loading  is  sometimes  called  the  fatigue  limit.  Cyclic  fatigue  in 
structural  materials  such  as  metals,  concrete  and  soil  has  been 
thoroughly  investigated  in  the  last  decades.  Rock  fatigue,  however, 
has  received  only  little  attention.  An  early  attempt  to  establish 
fatigue  characteristics  of  a  limestone  was  inconclusive  mainly  due  to 
a  very  limited  number  of  tests  undertaken  (3).  Burdine  (1),  however, 
ran  an  extensive  series  of  compressive  cyclic  loading  tests  in  Berea 
Sandstone,  discovered  that  the  rock  was  definitely  weakened  by  repet¬ 
itive  loading,  and  determined  S-N  curves  for  different  testing  condi¬ 
tions.  Hardy  and  Chugh  (4)  used  modern  equipment  to  improve  testing 
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thoroughly  investigated  in  the  last  decades.  Rock  fatigue,  however, 
has  received  only  little  attention.  An  early  attempt  to  establish 
fatigue  characteristics  of  a  limestone  was  inconclusive  mainly  due  to 
a  very  limited  number  of  tests  undertaken  (3).  Burdinc  (1),  however, 
ran  an  extensive  series  of  compressive  cyclic  loading  tests  in  Bcrc-a 
Sandstone,  discovered  that  the  rock  was  definitely  weakened  by  repet¬ 
itive  loading,  and  determined  S-N  curves  for  different  testing  condi¬ 
tions.  Hardy  and  Chugh  (4)  used  modern  equipment  to  improve  testing 
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procedures  and  concluded  that  three  additional  rock  types  were  fatigue 
prone  even  under  a  low  maximum  number  of  cycles  per  test  (10,  000). 

The  present  paper  reports  the  results  of  the  first  phase  in  a  compre¬ 
hensive  experimental  study  of  rock  cyclic  fatigue.  The  main  objectives 
of  the  investigation  are  to  provide  engineers  quantitative  results  in  the 
form  of  S-N  curves  obtained  for  different  rock  types  under  various 
stress  conditions,  to  determine  the  fatigue  strength  of  failed  rock  for 
the  use  of  mine  designers  and  earthquake  researchers,  and  to  provide 
a  better  understanding  of  the  fatigue  mechanism  by  strain  measure¬ 
ments  and  careful  observations  of  fabric  changes. 

LABORATORY  EQUIPMENT  AND  TEST  PROCEDURES 

Rock  Specimens 

Two  rock  types  have  been  used  in  the  reported  study.  White 
Tennessee  Marble  was  chosen  because  of  its  uniformity,  isotropy 
and  fine  grain  size.  It  has  been  extensively  tested  and  is  known  to 
have  very  consistent  mechanical  behavior  (6).  Georgia  Marble  was 
selected  for  the  ease  with  which  it  can  be  controlled  in  its  post  failure 
mode  (the  descending  part  of  the  complete  stress-strain  curve). 

Cylindrical  specimens,  1.  0  inch  in  diameter  and  2.  5  inch  long,  were 
cored  out  of  a  large  rock  block  by  diamond  drilling  in  one 
direction  only.  The  automatically-fed  coring  yielded  straight  smooth 
rock  cylinders  that  did  not  require  further  machining.  Specimen  ends 
were  surface  ground  until  flat  and  parallel  faces  were  obtained  to  with¬ 
in  0.  001  inch.  The  specimens  were  then  oven  dried  at  120°F  for  a 
week  prior  to  testing. 

Apparatus 


Specimens  were  loaded  in  an  electro-hydraulic  servo- controlled 


loading  machine  of  100,  000  lbs.  capacity.  A  general  view  of  the 
entire  apparatus  is  shown  in  F'ig.  1.  The  machine  can  be  programmed 


FIG.  1  LOADING  SYSTEM,  ELECTRONIC  CONTROLS  AND 
DATA  ACQUISITION  INSTRUMENTS. 

through  a  function  generator  to  apply  cyclic  loading  to  a  specimen  by 
controlling  either  stress  or  strain  rates.  In  the  experiments 
described  a  constant  rate  of  compressive  load  was  applied  to  speci¬ 
mens  using  a  triangular  wave  shape.  A  digital  counter  gave  the 
number  of  cycles  per  test,  and  an  X-Y  recorder  was  used  to  plot 
stress-strain  or  strain-time  curves.  The  load  applied  to  specimens 
was  sensed  by  a  dynamic  load  cell  mounted  on  the  bottom  part  of  the 
machine  crosshead.  The  longitudinal  strain  was  measured  through 
strain  gages  mounted  on  a  system  of  double  cantilevers  (Fig.  2). 

The  rings  holding  the  cantilever  device  were  threaded  on  to  the 
platens  in  contact  with  the  rock  specimens  and  held  in  place  by 
locking  nuts.  This  method  of  attachment  proved  superior  to  the  pre¬ 
viously  used  set-screws  which  tended  to  slip  in  long  duration  tests, 

3  A 


NOT  REPRODUCIBLE 


FIG.  2  LOADING  JIG  WITH  SPECIMEN,  CANTILEVER  SET 
AND  SWIVEL  HEAD. 

The  platens  had  a  diameter  only  slightly  larger  than  that  of  the 
specimens  and  were  part  of  a  loading  jig  especially  built  fot  these 
tests.  The  lower  platen  was  rigidly  attached  to  the  hydraulic  ram; 
the  upper  platen  contained  a  swivel  head  mechanism  to  ensure  that 
complete  contact  was  made  with  the  specimen.  To  prevent  slippage 
of  the  swivel  head,  six  bolts  were  used  to  lock  it  to  the  part  of  the 
loading  jig  rigidly  attached  to  the  crosshead.  The  locking  of  the  bolts 
was  a  tedious  job,  carefully  performed  to  ensure  good  alignment 
between  specimen  and  platens.  Overtightening  of  one  of  the  bolts 
could  cause  misalignment  and  hence  premature  failure  in  cyclic 
testing.  Utmost  care  was,  however,  taken  in  specimen  installation 


and  the  consistency  of  the  results  testify  to  the  success  of  the  method. 


Experimental  Program 

The  portion  of  the  fatigue  testing  program  that  has  been  carried 
out  so  far,  and  is  described  in  this  paper,  was  limited  to  uniaxial 
compression.  The  intent  was  to  test  rock  cylinders  under  com¬ 
pressive  loading  increased  steadily  from  zero  load  to  a  certain  upper 
peak,  decreased  at  the  same  rate  to  the  initial  point,  and  cycled 
until  failure  occurcd.  To  prevent  loss  of  contact  between  sample  and 
platens  at  zero  load,  the  lower  peak  of  the  compressive  cycle  was 
moved  up  to  around  200  psi  and  was  kept  approximately  constant 
throughout  the  testing  program.  The  tipper  peak,  however,  was  set 
anew  for  each  test. 

The  cyclic  rate  used  was  intended  to  simulate  the  frequencies  of 
the  major  pulses  in  earthquakes  (1-2  cps)  and  blasting  (  10  cps).  In 
the  tests  described  the  frequencies  were  kept  at  1-4  cps.  The  higher 
frequencies  were  preferred  because  of  shorter  test  duration,  but  the 
lower  frequencies  (1-2  cps)  were  necessary  when  stress- strain 
recordings  were  made,  due  to  the  X-Y  plotter  response  limitations. 

g 

The  maximum  number  of  cycles  per  test  was  kept  at  10  with  a  few 

exceptions  when  the  figure  was  exceeded.  To  obtain  a  representative 

S-N  curve  the  fatigue  life  of  specimens  (number  of  cycles  to  failure) 

was  determined  for  different  upper- peak  stress  values.  The  general 

procedure  was  to  reduce  the  maximum  applied  load  from  test  to  test 

6 

until  the  fatigue  life  of  the  specimens  reached  10  cycles. 

Cyclic  tests  were  also  run  on  specimens  that  had  been  first 
loaded  up  to  their  peak  compressive  stress  carrying  capacity.  In 
Georgia  Marble,  specimens  were  actually  brought  to  70%  of  the  peak 
stress  capacity  on  the  descending  side  of  the  complete  stress- strain 
curve,  and  only  then  cyclicly  loaded.  In  order  to  test  the  fatigue 
characteristics  of  these  failed  specimens  quasi-static  strain- 
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controlled  loading  was  used  to  bring  the  rock  to  the  initial  condition 
for  cyclic  stress- controlled  loading. 


EXPERIMENTAL  RESULTS 

The  most  exlcnsivc  testing  in  the  program  has  been  carried  out 
in  White  Tennessee  Marble.  The  first  stage  of  the  experimental  work 
was  to  determine  the  compressive  strength  of  the  rock,  a  value  to  be 
used  as  the  upper  limit  to  the  maximum  applied  stress  in  cyclic 
loading.  Specimens  were  loaded  stress- controlled,  at  four  different 
rates,  two  in  the  quasi -static  range,  and  two  in  the  dynamic  range 
used  in  the  cyclic  loading.  The  results  (Table  1)  show  that  the 

TABLE  1.-- COMPRESSIVE  STRENGTHS  OF  WHITE  TENNESSEE 
MARBLE  AT  DIFFERENT  LOADING  RATES 


Loading  rate 
(psi/sce) 

No.  of 
Specimens 

.  "  . . .  . .  " 

Mean  Compressive 

strength  (psi) 

Standard  Deviation 

j  psi 

% 

40 

Q 

o 

19, 740 

355 

1.  8 

100 

11 

21,150 

900 

4.  2 

50, 000 

12 

23,  285 

840 

3.  6 

200,  000 

7 

24, 890 

425 

1.  7 

compressive  sti  ength  varies  considerably  with  the  rate  of  loading. 
Hence,  whenever  this  parameter  is  employed  it  should  be  accompanied 
by  the  conditions  under  which  it  was  determined.  Because  of  the  high 
response  required, a  storage  type  oscilloscope  was  employed  to 
directly  record  the  stress  strain  curve  in  each  of  the  tests.  By 
photographing  the  oscilloscope  trace  a  permanent  record  was  obtained. 

The  first  important  conclusion  drawn  from  cyclic  testing  of  White 
Tennessee  Marble  is  that  it  is  definitely  weakened  by  repetitive  loading. 
The  fatigue  effect  can  best  be  verified  from  the  S-N  curve  shown  in 
Fig.  3.  The  stress  is  given  in  percentage  of  the  compressive  strength 

G/| 


at  50,  000  psi/sec.  Most  of  the  tests  were  run  at  4  cps,  while  some 
were  run  at  1  or  2  cps  for  the  purpose  of  recording  stress-strain 
characteristics.  No  apparent  difference  was  found  between  the 
frequencies  used  as  far  as  specimen  fatigue  life  was  concerned.  The 
S-N  curve  clearly  shows  that  as  the  maximum  compression  decreases 
the  life  expectancy  of  a  specimen  increases.  Fig.  3  also  indicates 
that  the  spread  of  experimental  points  was  surprisingly  limited,  due 
to  the  uniformity  of  the  rock  as  well  as  the  great  care  taken  in  the 
preparation  and  operation  of  the  tests.  In  a  semi-logarithmic  plot, 
as  the  one  shown  in  Fig.  3,  the  average  relationship  between  the 
maximum  applied  stress  (S)  and  the  number  of  cycles  needed  to  cause 
failure  (N)  is  given  by  a  straight  line.  At  75%  of  the  dynamic 
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compressive  strength  used  specimens  did  not  fail  within  the  experi- 

g 

mental  limit  of  10  cycles.  Moreover,  when  loaded  monotonically  to 
failure  at  the  conclusion  of  the  cyclic  test,  they  did  not  show  any  sign 
of  weakening  or  strengthening  effects.  It  is  hardly  expected  that  rock 
will  be  subjected  to  a  larger  number  of  similar  loading  cycles  during 
the  expected  life  of  an  engineering  structure.  In  an  earthquake,  for 
example,  no  more  than  100-200  cycles  are  encountered  (2).  Hence, 
the  value  of  17,450  psi  (-0.  75  x  23,  235  psi)  can  be  used  as  the  fatigue 
limit  or  the  critical  compressive  strength,  of  White  Tennessee  Marble. 
A  design  in  intact  rock  using  this  value,  not  only  is  protected  against 
static  and  dynamic,  stresses  but  also  against  all  kinds  of  cyclic  loading. 

Typical  stress  strain  cruvcs  are  shown  in  Fig.  4.  Permanent 


N  =  4  cycles  N=33  cycles  N  =  3700  cycles 


FIG. 


4  STRESS-STRAIN 
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strain  appears  to  accumulate  in  the  rock  with  the  dissipation  of  a 
relatively  high  amount  of  energy.  Common  to  all  tests  is  the  phenom- 
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enon  of  rather  large  hysteresis  in  the  first  few  cycles,  an  almost 
complete  closure  in  the  next  group  of  cycles,  and  a  reopening  in  the 
last  several  cycles  prior  to  failure.  The  hysteresis  in  the  first 
cycles  is  probably  due  to  the  very  high  peak  loads  which  take  the  rock 
to  beyond  its  linear  clastic  limit  and  into  the  zone  of  irreversible 
structure  changes.  In  the  very  short  life  tests,  the  damage  done  by 
the  high  stresses  is  irreparable  and  the  specimens  break  before 
reaching  the  second  stage.  For  lower  applied  maximum  stresses  the 
phenomenon  of  hysteresis  closure  occurs  seemingly  due  to  micro- 
crack  propagation  being  blocked  and  the  linear  elastic  limit  being 
lifted.  As  cyclic  loading  continues,  the  actual  fatigue  phenomenon 
eventually  takes  place,  internal  cracking  reinitiates,  the  hysteresis 
grows  and  the  specimen  finally  fails. 

The  amount  of  strain  difference  between  the  upper  peaks  of  the 
last  and  first  stress- strain  cycles  has  been  closely  followed.  No 
hard  conclusions  can  he  drawn  nr.  yet,  but  it  appears  that  it  is  defi¬ 
nitely  limited  by  the  strain  difference  between  the  ascending  and  the 
descending  parts  of  the  complete  stress-strain  curve  for  the  same 
value  of  stress.  This  implies  that  the  complete  stress-strain  curve 
actually  defines  the  limit  of  strain  that  can  he  applied  to  a  rock  at  a 
certain  stress  level,  without  producing  failure.  Quantitatively,  the 
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amount  of  strain  difference  varied  between  25  -  35  x  10  inch/inch, 
which  is  within  the  limit  set  by  the  complete  stress-strain  curve. 

The  behavior  of  strain  versus  time  was  often  recorded  and  a 
typical  result  is  shown  in  Fig.  5.  The  curve  shaped  by  the  upper  peak 
points  clearly  resembles  that  of  creep  behavior,  and  can  be  divided 
into  three  stages.  There  is  a  primary  stage  in  which  the  upper  peak 
strain  increases  at  a  decelerating  rale.  It  is  followed  by  a  steady 
state  stage  which  appears  in  Fig.  5  as  an  ascending  straight  line. 

In  medium  and  long  life  tests  this  stage  is  invariably  the*  longest.  The 
third  and  final  stage  is  that  of  accelerating  upper  peak  strain  culmi¬ 
nating  in  specimen  failure.  The  stages  observed  in  the  strain-time 


strain  characteristics. 


behavior  match  those  in  the  s  tress - 


FIG.  5  TYPICAL  STRAIN -TIME  CURVE 

It  can  be  easily  verified  from  Fig.  5  that  the ‘lower  peak  strain 

points  undergo  a  much  more  reduced  amount  of  increase.  This 

implies  that  the  value  of  the  tangent  modulus  per  cycle  generally 

decreases  during  the  test.  Indeed,  the  average  modulus  of  the 
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ascending  portion  of  the  first  cycle  was  11.  1  x  10  psi,  while  that  of  the 

Q 

last  cycle  was  9.  9  x  10  psi.  Similar  results  had  boon  found  in 
concrete  (5). 

Typical  fatigue  failed  specimens  arc  shown  in  Fig.  6.  No  apparent 
external  difference  was  found  between  this  type  of  failure  and  that 
encountered  in  quasi -static  loading.  Several  samples  wore  removed 
from  the  loading  machine  during  different  stages  of  the  cyclic  loading. 
The  samples  were  vertically  sectioned;  polished  and  photographed. 

No  damage  was  observed  in  specimens  that  had  not  reached  the  accel¬ 
erated  strain  increase  stage.  Those  that  had  reached  the  last  stage 
prior  lc)  failure  showed  extensive  structural  damage  dominated  by 
vertical  cracking.  A  photomicrograph  of  the;  central  part  of  a  vertical 
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section  in  a  specimen  removed  during  thelast  stage  is  presented  in 

Fig.  7. 


FIG.  6  FATIGUE  FAILED  WHITE  TENNESSEE  MARBLE  SPECIMENS 


A  series  of  cyclic  loading  tests  was  run  in  White  Tennessee  Marble 
that  had  first  been  loaded  in  strain  control  to  its  compressive  stress 
carrying  capacity  limit.  This  is  the  value  commonly  referred  to  as  the 
compressive  strength.  In  stress-controlled  loading  rock  specimens  will 
fail  violently  at  this  point  as  they  cannot  take  any  additional  compress  i\  o 
force.  However,  if  the  loading  is  performed  by  controlling  strain  rate 
the  rock  will  not  collapse  as  it  can  usually  continue  to  shrink  in  length 
while  its  load  supporting  capacity  is  gradually  lowered.  In  litis  fashion 
a  complete  stress- strain  curve  can  be  obtained.  Fig.  8  shows  such  a 
curve  in  addition  to  four  other  plots.  In  each  of  the  plots  a  specimen 
was  loaded  up  to  its  compressive  strength,  unloaded,  and  then  cycled 
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FIG.  7  PHOTOMICROGRAPH  OF  A  FATIGUE  FAILED  SPECIMEN- 

VERTICAL  SECTION 


in  stress  control  to  on  upper  peak  value  lower  than  the  ultimate 
strength.  The  purpose  was  to  verify  whether  such  failed  rock  could 
support  fatigue  type  loading.  Surprisingly  the  rock  appeared  rather 
strong,  and  could  probably  still  perform  useful  work,  although  a 
comparison  with  the  S~N  curve  in  Fig.  3  will  show  that  it  had  consi¬ 
derably  weakened.  As  expected,  the  number  of  cycles  increased  as 
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the  maximum  applied  load  decreased.  A  close  look  at  Fig.  8  will 
show  that  in  three  of  the  four  plots  the  upper  peak  of  the  last  cycle 
was  in  the  close  proximity  of  the  expected  descending  line  of  the 
complete  stress-strain  curve.  This  was  in  accordance  with  the 
hypothesis  raised  above  regarding  the  extent  of  the  upper  peak  strain 
increase  during  cyclic  loading. 

An  additional  series  of  tests  was  contemplated  well  within  the 
failed  region  characterized  by  the  descending  part  of  the  complete 
stress-strain  curve.  Due  to  the  very  steep  descending  stress-strain 
curve  in  White  Tennessee  Marble,  a  different  rock,  namely  Georgia 
Marble,  was  chosen  for  these  tests.  As  seen  in  Fig.  9  the  slope  of 
the  failed  portion  of  the  curve  is  very  mild  in  the  latter  rock.  This 
enabled  one  not  only  to  obtain  repeatable  complete  stress-strain 
curves,  but  also  to  stop  the  loading  at  any  point  in  the  failed  zone  and 
unload  without  the  danger  of  specimen  collapse.  Georgia  Marble 
specimens  were  loaded  in  strain  control  to  their  compressive  strength 
and  beyond, to  70%  of  the  load  carrying  capacity  within  the  failed  zone. 
They  were  then  unloaded,  and  cyclicly  loaded  in  stress  control  to 
upper  peak  values  of  different  magnitudes.  The  results  arc  shown  in 
Fig.  9.  A  preliminary  S-N  curve  based  on  these  tests  is  shown  in 
Fig.  10.  Again,  one  is  surprised  to  note  that  failed  rock  can  still 
show  remarkable  strength  as  far  as  cyclic  loading  carrying  capacity. 

It  is  emphasized  again,  however,  that  for  the  same  peak  loads  in  the 
unfailed  mode  the  number  of  cycles  would  be  appreciably  larger. 
Moreover,  it  is  clear  that  in  both  rocks. the  fatigue  limit  is  lowered  in 
failed  rock.  It  can  be  noted  from  Fig.  9  that  the  upper  peak  of  the  last 
cycle  falls  again  in  the  close  proximity  of  the  expected  descending 
complete  stress-strain  curve.  The  behavior  of  failed  rock  under 
cyclic  loading  is  of  particular*  importance  to  the  understanding  of 
failed  underground  pillars  and  walls  subjected  to  fatigue  type  stresses. 
Because  failed  rock  is  internally  fractured,  the  study  may  also  be 
indicative  of  the  behavior  of  jointed  benches  or  slopes  and  faulted 


formations  subjected  to  earthquake,  blasting,  or  traffic  loadings. 


FIG.  9  STRESS-STRAIN  CURVES  FOR  FAILED  GEORGIA  MARBLE 

(DASHED  LINES  ARE  EXPECTED 
COMPLETE  STRESS-STRAIN  CURVES) 


NUMBER  OF  CYCLES  (N) 

FIG.  10  S-N  CURVE  FOR  FAILED  GEORGIA  MARBLE 

SUMMARY  AND  CONCLUSIONS 

The  work  reported  here  is  only  the  first  phase  of  an  extensive 
investigation  into  the  behavior  of  rock  under  cyclic  loading.  The  main 
results  obtained  finis  far  in  cyclic  uniaxial  compression  of  White 
Tennessee  Marble  are: 

1)  Cyclic  loading  has  a  definite  weakening  (fatigue)  effect  on  the 
rock  if  the  maximum  applied  stress  is  in  the  75-100%  range  of  the 
compressive  strength. 

2)  Stress- strain  records  reveal  three  major  stages  occuring 
during  cyclic  loading:  decreasing  hysteresis,  no  hysteresis,  increasing 
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hysteresis.  The  Inst  stage  is  associated  with  fatigue  affected  internal 
fracturing. 

3)  The  increase  in  upper  peak  strain  between  the  first  and  the 
last  cycles  is  always  within  the  limit  set  by  the  complete  stress- strain 
curve  for  the  same  maximum  applied  load. 

4)  The  variation  of  the  cyclic  upper  peak  strain  with  time  closely 
resembles  that  of  creep,  with  primary,  steady  state,  and  tertiary  or 
accelerated  stages. 

5)  No  visible  structural  damage  can  be  observed  in  specimens 
removed  prior  to  the  accelerated  strain-increase  stage.  However, 
polished  sections  of  specimens  that  had  reached  the  tertiary  stage 
show  an  abundance  of  internal  fracturing  dominated  by  vertical 
cracking  along  and  across  grain  boundaries. 

G)  Failed  rock,  although  weakened,  can  sustain  a  certain  amount 
of  fatigue  loading  depending  on  the  level  of  applied  stresses.  This 
conclusion  is  supported  by  testing  of  both  Tennessee  and  Georgia 
Marbles. 

The  results  obtained  thus  far  show  that  rock  cyclic  fatigue  is 
a  phenomenon  that  cannot  be  ignored  by  engineers  and  rock  mechanics 
scientists.  A  direct  recommendation  emerging  from  the  study  is  that 
in  surface  and  underground  design  the  apparent  fatigue  limit,  as 
determined  by  lests  similar  to  those  described  above,  be  used  instead 
of  the  commonly  employed  compressive  strength  value.  In  the  design 
of  structures  that  eventually  reach  their  peak  load  carrying  capacity, 
use  could  be  made  of  (he  fatigue  characteristics  of  failed  rock.  These 
characteristics  arc  believed  to  also  be  indicative  of  jointed  and  faulted 
rock  behavior  under  cyclic  loading.  The  mechanism  of  fatigue  is  not 
yet  clearly  understood  but  the  findings  arc  promising  enough  to  justify 
continuation  of  the  research.  Studies  on  the  behavior  of  rock  in  uni¬ 
axial  compression,  tension  and  triaxial  compression  subjected  to  cyclic 
fatigue  arc  currently  underway. 


17  // 


A  CKN'OW  Ll'lDGK  i\l  RNTS 


The  authors  wish  (o  acknowledge  (he:  financial  assistance  of  the 
Wisconsin  Alumni  Research  Foundation  in  the  early  stages  of  the 
investigation.  The  research  was  supported  by  the  Advanced  Research 
Projects  Agency  of  the  Department  of  Defense  and  was  monitored  by 
the  United  Stales  llurcau  of  Alines,  Twin  Cities  Re;  search  Center, 
under  contract  no.  I  lOHlOOOv! .  .1.  A.  Hudson  and  F.  Rummcl, 

University  of  TU innesola,  provided  invaluable  suggestions  and  assistance. 


APPENDIX-- RFFERFNCFS 


1.  Burdine,  ]s’.  T. ,  "Rock  Failure  l'mlcf  Dynamic  Loading 
Conditions",  Society  of  Petroleum  En  mjaeoru  Journal,  Alar  eh 
19G3,  pp  1-8. 

2.  Ellis,  Willard  and  Hartman,  V.  T3. ,  "Dynamic  Soil  Strength  and 
Slope  Stability",  Journal  of  the  Soil  Mechanics  and  Foundation?; 
Division,  A.  S.  C.  E.  ,  vol.  83,  No.  SiVb],  July  1087,  pp.  355-37D. 

3.  Grover,  II.  J. ,  Dchlinger,  P.  ,  and  McClure,  G.  M.,  "Investi¬ 
gation  of  Fatigue  Characteristics  of  Hocks",  Report  by  Baltclle 
Memorial  Institute  to  Drilling  Research,  Incorporated,  Nov.  30, 
1980. 

4.  Hardy,  11.  H.  ,  Jr.  and  Chugh,  Y.  P. ,  "Failure  of  Geological 
Materials  Under  Low- Cycle  Fatigue",  Proceedings  of  the  Sixth 
Canadian  Symposium  on  Rock  Mechanics,  Montreal,  Canada, 

May  1970. 

5.  Nortlby,  Gene  M.  "Fatigue  of  Concrete- -A  Review  of  Research", 
Journal  of  the  American  Concrete  Institute,  Aug.  1958,  pp.  191-219. 

G.  Wawersik,  Wolfgang  R. ,  "Detailed  Analysis  of  Rock  Failure  in 
Laboratory  Compression  Tests",  thesis  presented  to  the 
University  of  Minnesota.  Minni-apolis,  in  19GG,  in  partial  fulfill¬ 
ment  of  the  requirements  for  the  degree  of  Doctor  of  Philosophy. 


